Electrophoresis of DNA on a disordered two-dimensional substrate 
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We propose a new method for electrophoretic separation of DNA in which adsorbed polymers are 
driven over a disordered two-dimensional substrate which contains attractive sites for the polymers. 
Using simulations of a model for long polymer chains, we show that the mobility increases with 
polymer length, in contrast to gel electrophoresis techniques, and that separation can be achieved for 
a range of length scales. We demonstrate that the separation mechanism relies on steric interactions 
between polymer segments, which prevent substrate disorder sites from trapping more than one 
DNA segment each. Since thermal activation does not play a significant role in determining the 
polymer mobility, band broadening due to diffusion can be avoided in our separation method. 

PACS numbers: 87.15.Tt, 87.14.Ge, 87.15.Aa 
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I. INTRODUCTION 

The development of new methods for efficiently sepa- 
rating charged biopolymers by length has been an area of 
significant recent activity due to the fact that strategies 
for genome sequencing are based on sorting DNA frag- 
ments by size |l| . Simple charge-based sorting is not pos- 
sible because the increase in electrostatic force on longer 
molecules with greater total charge is exactly offset by 
a corresponding increase in hydrodynamic drag Q. In- 
stead, separation is achieved using techniques such as 
gel electrophoresis, in which longer molecules are slowed 
relative to shorter ones due to interactions with cross 
links in the gel. Gel and capillary electrophoretic tech- 
niques are limited to DNA strands less than 4 x 10 4 base 
pairs (bp) in length §; the mobility saturates for longer 
strands, and sufficiently large strands fail to pass through 
the gel at all. There is a need for separation of strands 
up to 1 x 10 6 bp, and thus new techniques which can 
sort longer molecules are of particular interest Q. The 
motion of elastic strings through random media is also 
of general interest for a wide range of systems including 
magnetic domain wall motion, vortex lattice motion in 
superconductors, and charge density waves. 

Several recent proposals for electrophoretic techniques 
move away from the traditional media of gels and poly- 
mers and instead take advantage of advances in nano- 
lithographyto create microstructured devices for separa- 
tion g Ell HSHGIlHEEim. The sorting 

effectiveness of these techniques is limited by the relative 
size of the nanofabricated structure and the polymers to 
be sorted, making it necessary to fabricate a separate de- 
vice for each size range of interest. In contrast, Seo et 
al. 0| proposed an adsorption-based separation tech- 
nique that could permit the sorting of polymers which 
vary in size by three orders of magnitude. When the 
ionic strength of the buffer solution is altered the 
DNA is partially adsorbed onto a clean surface, form- 
ing a series of loops which extend into the solution and 
trains which are adsorbed on the substrate. It has been 
proposed that separation occurs because the longer poly- 



mers have a larger number of train segments, and thus 
experience a greater retardation of their motion [Tsl Il9| . 

The quasi-two dimensional geometry considered in 
Ref. |l6l | is very appealing for separation purposes, in 
part because adsorb ed p olymers spread out significantly 
into flat 'pancakes' |2(| compared to their coiled three- 
dimensional configurations, permitting better coupling to 
length differences. The sorting mechanism in Ref. |l6| 
precludes complete adsorption, however, since there is no 
separation for fully desorbed or fully adsorbed polymers. 
This limits the length range that can be processed, since 
if the surface is strongly attractive to DNA, long DNA 
chains fully adsorb and separation by length is lost. If 
instead the surface weakly attracts DNA, short chains 
desorb from the surface and cannot be separated pi) . 

Here we propose an alternative sorting technique for 
long DNA strands in which the polymers are fully ad- 
sorbed on the surface. To permit separation, we spatially 
modify the surface, but instead of of using posts or other 
impenetrable barriers, we consider randomly spaced pin- 
ning sites which temporarily retard the motion of the 
polymer, yet still allow it to pass through. Such pinning 
could be created via the manipulation of lipid bilayer 
membranes [22] or surface patterning |2^| . We show that 
in this geometry, longer polymers are more mobile than 
shorter ones, in contrast to typical separation methods 
where longer polymers move more slowly. This avoids 
the jamming or clogging associated with long polymers 
in other techniques. The steric interaction between poly- 
mer segments causes the longer polymers to be less well 
pinned by the random disorder than the short polymers, 
and allows separation by length to occur. 

To demonstrate our separation mechanism, we use a 
simulation model that we have developed for long DNA 
fragments. Many of the existing simulation models for 
electrophoretic processes are best suited for shorter poly- 
mers [24| . Since we are concerned with polymers up to 
300 /im in length, we do not attempt to simulate each 
atom in the polymer. Instead, we adopt a bead-spring 
model in which the polymer is represented by multi- 
ple beads |25j which are each spaced many persistence 
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lengths apart. There is an entropic resistance to the 
stretching of the polymer segment between two beads, 
which is represented by a finitely extensible nonlinear 
spring (FENE) potential [26j that replaces the internal 
degrees of freedom of the polymer molecule |27j. An es- 
sential assumption of this model is that the polymer seg- 
ment between beads is significantly longer than the poly- 
mer persistence length. This is in contrast to bead-stick 
models [28j, where the distance between beads is ten or 
less actual chemical segments. 



II. SIMULATION 

We employ Brownian dynamics |29|. permitting us to 
use time steps of order 0.1 ns, orders of magnitude greater 
than the sub-fs time steps required in all-atom molecular 
dynamics. In this technique, the solvent is treated sta- 
tistically rather than explicitly 30] . The dimcnsionlcss 
force on bead i in a chain L base pairs long represented 
by N beads is given by 
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where F FENE is the spring force along the chain, F EV 
represents the excluded volume between beads, F s is 
the force from a disordered substrate, F E is the elec- 
trophoretic force, and F T is a thermal noise term. Dis- 
tances are measured in terms of a s , the root mean square 
length of the spring. Forces are expressed in terms of 
ksTjas. We can neglect hydrodynamic interactions since 
they are screened due to the proximity to the solid sub- 
strate [jOI ■ Electroosmotic effects can be controlled 
in the usual way by means of a high concentration buffer 
[Hi El EJ HH . We assume that the Debye length is con- 
siderably smaller than the distance between the beads in 
our model. 

The force between bead i and neighboring beads is 
given by 
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where Q = 1 1 — lo | is the elongation of the spring, 1 
is the distance vector between two neighboring beads, 
lo = o"sl is the equilibrium spring length, Qo is the maxi- 
mum allowable elongation, and the Hookean spring con- 
stant H = 3/er 3 . This phenomenological spring potential 
[2f| has the properties that it is equivalent to a Hookean 
spring for small Q, but becomes infinite at finite spring 
elongation. The persistence length l p of double-stranded 
DNA is l p w 500 10. The Kuhn length b k = 2l p , giving 
frfc = 0.1/im, where we assume that the ionic strength of 
the buffer solution is sufficient to screen electrostatic re- 
pulsion between sections of the chain [^J . Each base pai r 
is 0.34 nm long so one Kuhn length contains 300 bp |3fil |. 
Since we will be using a Gaussian chain model for the 



excluded volume interactions, which requires the chain 
segments between beads to be represented statistically 
and not deterministically, the number n of Kuhn lengths 
between beads must be sufficiently large (well above the 
bead-rod limit of n — 1), so we take Qo = nbk = 1.6/im. 
This gives n — 16 and a s — y/nbk = 0.4/im. 

The excluded volume term accounts for the repulsive 
interaction between polymer segments when they ap- 
proach each other (37L l38| , an effect which is more pro- 
nounced in two dimensions than in three dimensions |39j . 
The excluded volume interaction for beads i and j a unit- 
less distance r b b apart is taken after that used in Ref. ^3 j 
which is based on the energy penalty due to overlap of 
two Gaussian coils, and has the form 



t?EV 



-Ar bb e Br >»>r b b 



where 



5/2 



3a 2 



(3) 



(4) 



Here the size parameter Si = nb\/Q 40], while the ex- 
cluded volume parameter v is taken to be v = 6|. This 
gives 



.4 



243V2^ 



B 



(5) 



The substrate roughness is represented by N p finite 
range parabolic pinning traps of radius a p = 0.4cr s , 
strength / p , and density p p . The pin size was chosen 
to be close to the Kuhn length. The force on bead i from 
pin k a distance r bp away is given by 



ik 



(6) 



where is the Heaviside step function. 

The electrophoretic force on each bead from an applied 
electric field E is 



F = qEy 



(7) 



Here, q — Xnb^ is the charge per bead, where the charge 
per unit length of DNA in solution is A = 4.6 x 10~ 10 
C/m, or O.3e^/A0. F T is the Langevin thermal noise 
term representing the Brownian forces. It is a delta- 
correlated white noise process which obeys (F T ) = and 
the fluctuation-dissipation theorem 41] {Fj \t)Fj \t + 
At)) = 2kBT(^ 1 5ij5(Ar). Here is a Kronecker delta 
tensor and <5(At) is the Dirac delta function. Time 
is measured in units of r = £<t 2 /(/cbT) and we take 
St = 0.001. C is the friction coefficient characterizing 
the viscous interaction between the bead and the solvent. 
We use the experimentally measured value of £ for poly- 
mers diffusing in a bilayer, £ = 2.97 x 10~ 7 Ns/m |42| . 
Theoretically, for a polymer moving in three dimensions, 
C = 67r?7 s <7, where a is the effective bead radius and r] s is 
the solvent viscosity. We note that a theoretical expres- 
sion for the friction coefficient in the case of a particle 
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FIG. 1: Average velocity (V) versus electric field for pinning 
density p p — 0.8, pin strength f p = 60, and polymers of 
different length represented by the number of beads N = 
5 (circle), 10 (square), 20 (diamond), 30 (triangle up), 40 
(triangle left), 50 (triangle right), 80 (triangle down), and 
100 (plus). 



confined to a membrane suspended in a solvent has been 
developed by Saffman |43j, where only a weak depen- 
dence of £ on effective bead radius is obtained. 



III. RESULTS 

We first consider the velocity of the polymers over the 
rough substrate as a function of polymer length L. We 
sweep the electric field strength and find the average ve- 
locity (V) at each field value during 200 repetitions of 
the sweep. In Fig. \T\we plot the velocity- force curves for 
polymers of length ranging from ./V = 5 to 100, where 
N is the number of beads used to model the polymer, 
in a sample with pinning density p p — 0.8 and strength 
f p = 60 at room temperature. In physical units, this 
length range is 8^m to 160/im, and it includes A-phage 
DNA, which has a contour length of 21.2/xm 0]. A ve- 
locity of 10 corresponds to 0.3 /im/s, and an applied field 
of 30 corresponds to 4.16 V/cm. After the polymers de- 
pin, there is a range of driving force over which we find a 
nonlinear velocity-force characteristic. Within this range 
the shorter polymers move more slowly than the longer 
polymers for a given electric field strength. 

Short polymers are better pinned by the underlying 
disorder than long polymers, and thus a higher driving 
force must be applied before the short polymers begin 
to move over the substrate. In Fig. [21 we illustrate the 
dependence of the critical depinning force F c on polymer 
length L for a range of pinning strengths and densities. 
Here we define F c as the driving force at which (V) = 1. 
We have chosen this definition since in our technique, 
physical separation of the polymers will occur when the 
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FIG. 2: (a) Critical depinning force F c versus polymer length 
L on a linear-log scale for fixed f p — 60 and pinning density 
p p = 0.2 (circles), 0.4 (squares), 0.6 (diamonds), 0.8 (triangle 
up), and 1.0 (triangle left). Dashed lines are logarithmic fits, 
(b) Scaled critical depinning force F c / f P versus L for fixed 
p p = 0.8 and f p = 20 (circles), 40 (squares), 60 (diamonds), 
and 80 (triangles). Dashed lines are logarithmic fits. 



polymers are moving with different velocities. In an ex- 
periment, if the electric field were held at a value between 
F c for polymers of two lengths, the two polymers will be 
separated since one is strongly mobile while the other is 
nearly immobile. As shown in Fig. [21 in each case we 
find that F c drops logarithmically with L, as indicated 
by the dashed lines. As the pinning density is reduced 
from p = 1.0 to p — 0.2, shown in Fig. Ufa), the depin- 
ning force drops and the variation of F c with L becomes 
steeper, meaning that the separation resolution is en- 
hanced. At the same time, the length range over which 
effective separation can be achieved drops, and thus there 
is a trade-off which must be considered depending on the 
range of sizes that are to be separated. We show the 
scaled length dependence of F c /f p in Fig. |2tb) for pin- 
ning strengths f p — 20 to 80. The pinning effectiveness 
drops slightly faster than the pinning strength, as indi- 
cated by the fact that the curves do not fall on top of each 
other. For the weakest pins, the separation effectiveness 
washes out above a length of N ~ 100. 

The excluded volume interactions play the key role in 
the separation mechanism. What is happening physically 
can be understood as follows. Consider a polymer com- 
posed of only a single bead. This polymer can be com- 
pletely trapped by a single pinning site. Next, consider a 
polymer composed of two beads. Although it is possible 
for the polymer to find two adjacent pinning sites such 
that both beads are pinned, it is more likely that one bead 
will be pinned while the other is still free. In this case the 
force from only one pin will have to hold two beads still 
against the electrophoretic force. If the excluded volume 
interactions are removed, both beads can fit inside the 
pin, doubling the effective pinning force. As the poly- 
mer becomes longer and is composed of more beads, the 
relative fraction of the polymer that is pinned decreases, 
provided that each pin can capture only one bead. This 
results in the decreased threshold for depinning and the 
increased mobility of the longer polymers relative to the 
short ones. The importance of the excluded volume in- 
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FIG. 3: Average velocity (V) versus electric field for a sys- 
tem with no excluded volume at pinning density p p = 0.8, pin 
strength f p — 60, and polymers of different length N = 5 (cir- 
cle), 10 (square), 20 (diamond), 30 (triangle up), 40 (triangle 
left), 50 (triangle right), 80 (triangle down), and 100 (plus). 
Inset: F c versus L for the same system. 



teraction is that it enforces a pin occupancy of at most 
one bead per pin. Thus, the excluded volume interaction 
is what produces the decrease in F c with polymer length. 
We test this by running a series of simulations without 
the excluded volume interaction. The depinning force F c 
for this case is shown in the inset to Fig. |3 where it is 
clear that F c has no significant dependence on L. The 
corresponding velocity-force curves are shown in Fig. |3 
where it can clearly be seen that the curves lie on top of 
each other except for the very shortest polymers. 

We stress that the separation mechanism at work here 
is significantly different than that which occurs in the 
case of impenetrable obstacles such as cross links in gels 
or nanofabricated posts. This can be seen by observing 
images of the moving polymers. A representative set of 
images for polymers of different length is shown in Fig. ^ 
The chains are moving toward the top of the figure in the 
+y direction. Rather than forming hairpin structures, 
the polymers frequently form a bundle on their advancing 
end, and sometimes drag one or two tail segments. 

For separation purposes, the polymer velocity must de- 
pend on length. This can be achieved if the depinning 
force of the polymers is length dependent. To demon- 
strate this explicitly, we run a series of simulations in 
which the driving force is held at a fixed value, and mea- 
sure the average velocity (V). The results are plotted in 
Fig. for five different values of F E in a sample with 
p p = 0.8 and f p = 60. The velocity increases logarith- 
mically with polymer length, as indicated by the dashed 
lines, and velocity variations of an order of magnitude 
can be achieved. 

We note that in traditional gel electrophoresis tech- 
niques, diffusion is an important limiting effect, since the 




(f) (g) (h) (i) (j) 

FIG. 4: Images from polymer simulation showing bead posi- 
tions for length N = (a) 40, (b) 50, (c) 60, (d) 70, (e) 80, (f) 
90, (g) 100, (h) 125, (i) 150, (j) 175, and (k) 200. The driving 
force is in the +y direction, toward the top of the figure, as 
indicated by the arrow in (a). 
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FIG. 5: (a) Velocity (V) versus polymer length L for f p = 60 
and p p = 0.8 at different applied driving fields of F E = 37.88 
(circles), 39.0 (squares), 40.12 (diamonds), 41.24 (triangle 
up), and 42.36 (triangle left). Dashed lines are logarithmic 
fits. 
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FIG. 6: (a) Distance d in units of a s traveled over a fixed time 
interval versus polymer length L for f p = 60 and p p = 0.8 at 
F E = 37.88 (circles), 39.0 (squares), 40.12 (diamonds), 41.24 
(triangle up), and 42.36 (triangle left). Error bars indicate 
the spread in distance traveled over 100 realizations of disor- 
der, (b) Distance d traveled over a longer fixed time interval 
versus polymer length for the same system as in panel (a) 
with F E = 37.88 (circles), 39.0 (squares), 40.12 (diamonds), 
41.24 (triangle up), and 42.36 (triangle left). Error bars in- 
dicate the spread in distance traveled over 20 realizations of 
disorder. 



polymers are moving through the gel relatively slowly 
and depend on thermal fluctuations to help them translo- 
cate through the gel. This effect is particularly pro- 
nounced for long polymers, which have extreme difficulty 
passing through the gel at all. In contrast, in the tech- 
nique proposed here, the polymers are much more mo- 
bile than they would be in a gel. The configurations and 
depinning of the polymers are dominated by the strong 
electric fields and pinning imposed, and thermal effects 
play essentially no role in the separation. We observe 
no significant thermal diffusion in our system at all. As 
a result, diffusive broadening of the bands can be pre- 
vented. The bands do still broaden due to the intrinsic 
randomness of the pinning, which causes the progress of 
the polymers over the substrate to be somewhat vari- 
able. To illustrate the magnitude of this broadening, in 
Fig. EJa) we plot the total distance traveled by the poly- 
mers under different drives applied for a fixed period of 
time for 100 realizations of disorder. Error bars indicate 
the average maximum and minimum distances traveled 
by polymers of a particular length. Higher resolution can 
be obtained by allowing the polymers to move a larger 
distance through the gel, as illustrated in Fig. EJb). 

We quantify the separation power of our technique by 
measuring the resolution as a function of polymer length. 
The resolution is affected by both the selectivity and effi- 
ciency of the separation for a given length difference . 
The selectivity AV/V m is proportional to the difference 
in mobility for polymers of different length, 



AVj = 2{{V l+1 ) - (Vj)) 
V m {V i+ i) + (K) 



(8) 



where (V*) is the average velocity for a polymer of length 
Li- Fig. [7{ a) shows that the selectivity for the same 
system in Fig.|SJa) does not vary with F E and is highest 
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FIG. 7: (a) Selectivity AV/V m from the system in Fig. EI a) 
versus polymer length L for f p — 60 and p v — 0.8 at F E = 
37.88 (circles), 39.0 (squares), 40.12 (diamonds), 41.24 (tri- 
angle up), and 42.36 (triangle left), (b) Efficiency Ne for 
the same system at F E — 37.88 (open circles), 39.0 (filled 
squares), 40.12 (diamonds), 41.24 (filled triangle up), and 
42.36 (triangle left), (c) Resolution R for the same system 
at F E = 37.88 (open circles), 39.0 (filled squares), 40.12 (di- 
amonds), 41.24 (filled triangle up), and 42.36 (triangle left), 
(d) Resolution at F E — 40.12 for: (open diamonds) the sys- 
tem in Fig. EI a) and (filled diamonds) the system in Fig.EJb) 
where the polymers were allowed to travel a longer distance. 



for the shortest polymers, in the same region where Fig.[S] 
indicates that the (V) versus L curve has the steepest 
slope. The efficiency Ne is proportional to the width of 
the band a^. observed after the polymers have traveled a 
distance x, 



(9) 



As can be seen in Fig. Efb), Ne increases with both L 
and F E , consistent with the decrease in the size of the 
error bars at higher L and F E shown in in Fig.UJa). The 
resolution R is defined as 



R 



(10) 



where Ne is the mean efficiency for the polymer lengths 
being compared. We plot R versus L in Fig. |7fc). The 
resolution depends more strongly on the selectivity than 
on the efficiency, and as a result we find that R is high- 
est for the shortest polymers and is not a strong function 
of F . The resolution can be improved by allowing the 
polymers to travel a longer distance, as in Fig.[|Jb). We 
compare the resolution for shorter and longer distances 
traveled in Fig.|J{d), where we find not only an enhance- 
ment of R for the longer travel distance, but also a shift 
in the peak value of R towards longer polymers. This 
suggests that the technique could be optimized for sepa- 
ration of the desired range of L by adjusting the distance 
traveled by the polymers. 



6 



IV. CONCLUSION 

In summary, we have used a model developed for the 
simulation of long DNA segments to demonstrate a new 
length separation mechanism for polymers adsorbed to a 
disordered two-dimensional substrate. Longer polymers 
are more mobile than short polymers, and the depin- 
ning force decreases logarithmically with polymer length. 
Correspondingly, the polymer velocity increases logarith- 
mically with length. The separation mechanism arises 
due to the excluded volume interaction between chain 
segments, which serves to reduce the effectiveness of the 
random pinning for longer polymers. One possible exper- 
imental system in which our proposed separation mech- 
anism could be realized is solid-supported cationic lipid 
membranes, where DNA is confined to two dimensions 
but free to diffuse in plane |3l| . The pinning could be pro- 
duced in the form of disorder on the supporting substrate, 
which would perturb the bilayer and interfere with the 
free diffusion of the DNA. Such disorder could potentially 
be produced by an experimental technique as simple as 
not fully cleaning the substrate before depositing the bi- 
layer. Our proposed separation mechanism offers several 



advantages over existing techniques. (1) It may not be 
necessary to use elaborate nanofabrication methods to 
produce the pinning. (2) The technique can be used to 
separate extremely long strands of DNA which will not 
pass through conventional gels. (3) It may be possible 
to achieve high throughput since the polymers do not 
need to work their way around fixed impassible obstacles, 
but are instead only temporarily hindered by the pinning 
sites, and can thus achieve much higher overall mobilities 
than are possible in a gel, particularly for long polymers. 
(4) Since thermal effects do not play a significant role 
in the separation technique, thermal broadening of the 
bands by diffusion should be strongly suppressed. The 
resolution limitation caused by band broadening from the 
intrinsic disorder of the substrate can be reduced by al- 
lowing the polymers to travel a longer distance during 
separation. 
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